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1 Introduction 



The function smx/x is endlessly fascinating. By setting x = tt/2 in the infinite product expansion 

oo 

sin a; n x , 

k=l 

one gets the first actual formula for ir that mankind ever discovered, dating from 1593 and due to 
Frangois Viete (1540-1603), whose Latinized name is Vieta. (Was any notice taken of the formula's 
400th anniversary, perhaps by the issue of a postage stamp?) From the samples of a function f(x) 
at equally spaced points x n , n £ Z, one can reconstruct the complete function with the aid of 
smx/x, provided / is "banddimited" and the spacing of the samples is small enough. This is the 
content of the Sampling Theorem, which lends its name to smx/x as the sampling function. 
Its importance in signal processing, where it is also known as sine x, is the result of its Fourier 
transform being the characteristic function of the interval [—1,1] (modulo a scalar factor). 



In section 2 we prove the infinite product expansion for sinx/x and derive Viete's formula. In 
section 3 we transform the product expansion with the Fourier transform and use convolution and 
delta distributions to prove it in a way that reveals a host of similar identities. Section 4 puts these 
identities into a probabilistic setting, and in section 5 we alter the probability experiments in order 
to make connections between infinite cosine products, Cantor sets, and sums of series with random 
signs, particularly the harmonic series. This leaves us with some interesting unsolved problems and 
conjectures for further work. 

"This article has appeared in Amer. Math. Monthly, 102:716-724, 1995. 
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2 An Elementary Proof 



Repeated use of the double angle formula for the sine shows that 



But 



X X 

sin x = 2 sin — cos — 
2 2 

XXX 

= 4 sin — cos — cos — 
4 4 2 



k=i 



lim 2 n sin — = x, 

n— +oo 2 n 



thereby proving the identity. See Figure 1 for an indication of how quickly the product converges. 
Let x = ir/2, make use of the half-angle identity, and there you have Viete's formula for n, 



2 V2 y/2 + V2 V 2 + V2 + y/2 



TT 



(2) 



At this point the cosine identity could remain an isolated curiosity of historical interest, relegated 
to the ends of exercise sets in textbooks. In fact, it is just the first of an infinite family of cosine 
product identities for s'mx/x. 




Figure 1. Graphs of s'mx/x and cos § cos | cos |. 
Where both graphs are visible, sinx/x is nearer the x axis. 



3 The Fourier Transform and More Identities 



For a complex valued function f(x) defined on the real line, the Fourier transform puts together / 
as a continuous linear combination of the "pure" oscillations e tuJX in which the coefficient in front 



2 



of e tuJX is denoted by f(oj). Thus, 

/oo 
/He^da;. (3) 
-oo 

The function / is the Fourier transform of / and the integral above is a description of how to 
get back / from / and is actually the formula for the inverse transform. How do we get / from /? 
That is given by this integral: 

/H = ^ J~ f(x)e-™dx. (4) 

Of course, the proofs of these relationships involve hypotheses on the functions so that the inte- 
grals make sense, but they can be extended beyond the realm of ordinary functions to generalized 
functions or distributions. We need more than ordinary functions in order to make sense of the 
Fourier transform of a sine or cosine. 



Notation: we also write the Fourier transform of / as J-(f) and the inverse transform of (ft as 

Consider cos&x, which by Euler's Identity may be written as 

cosbx = ^(e ibx + e~ ibx ). 

This shows the function written as a linear combination of just two of of the functions e tuJX for 
oj = b and u> = —b. The coefficients appear to be 1/2, but if we use them in the integral form 
with all other coefficients zero, then we cannot represent the cosine function. Instead, we must 
regard the coefficients as point masses at b and —b. Therefore, the Fourier transform of cos bx is 
(1/2) (Sb + S-b), where 5b denotes the Dirac delta distribution or point mass at the point b. All of 
this can be made rigorous, but at the expense of some long development in graduate level analysis. 
The approach here is at about the level of a second year course in engineering mathematics. 

In addition, the Fourier transform behaves nicely on a product of functions and turns it into the 
convolution of the transforms: 

fg = f* g- (5) 

For two functions (ft(oj) and ip(oj), the convolution (ft * tp is defined by 

/oo 
(ft(a)ift(uj — a) da. (6) 
-oo 

Again, we must extend convolution beyond the realm of functions. In particular we need convo- 
lutions of delta distributions and for them we can easily show that So behaves as the identity for 
convolution 

So * (ft = (ft (7) 

and that 

5 a * 5b = 5 a+b . (8) 



Now back to the cosine identity. Let f(x) = YlkLi cos(x/2 k ) and let f n be the nth partial product. 
The Fourier transform of /„ is 

" 1 

fn = * JJ 2^1/2" + <*-l/2")- 
fc=l 
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The asterisk in front of the product sign indicates a repeated convolution of the factors. Expanding 
for n = 3 we see that 

h = o (<*-7/8 + 6-5/8 H 6 7/8 ). 



Likewise 



2" z 4 ' 

beB„ 



where B n is the set of 2 n equally spaced numbers from — l + l/2 n to 1 — 1/2™ with spacing 
2/2™ = 1/2(™ _1 ). 

The sequence of measures /„ converges to the uniform density on [—1, 1] of total mass 1, which we 
can write as (l/2)x[-i,i]dw- The inverse transform is easy to compute: 



fj 



e^ x duj= Smr 



2 x 

The spectrum of (s'mx)/x is uniform in the interval — 1 < u> < 1. This means that sinx/x is a 
continuous linear combination of the "pure" harmonics e tu)X with the same weight of 1/2 for each 
w€[-l,l]. 

With this proof we have a way to generate a family of similar identities. Let us put point masses 
at 3 n equally spaced points from — 1 + 1/3™ to 1 — l/3 n with spacing 2/3™. Such a measure is the 
convolution *Hfc=i |(6-2/3 fc + 6o + 6 2 / 3 fc). Applying the inverse transform 

Q (<*-2/3* +6o + 6 2 /3fe)) = \ ( 2cos p + 1 
and taking limits gives us the infinite product identity 

(9) 



nl / 2x\ sin: 

n=l x 7 



Let us use the positive integer p as the base (we have just seen p = 2 and p = 3). The first measure 
/i is the sum of point masses at p points equally spaced from —1 + 1/p to 1 — 1/p with spacing 2/p. 

/i = -(tfi^ +53^ + <5^p + --- + <y £ -i) (10) 
p p p p p 

1 

= -V<w P (ii) 

y 1=0 

We let 



n / p-1 
fn = * TT - V] 6 2 i+i- P 



k=l \ P 1=0 

and one can see that f n consists of p n point masses equally spaced from —1 + l/p n to 1 — l/p n with 
spacing 2/p n . Taking the inverse transform we see that 

n p—1 



= n -E^^ + 1 -p) ix /p k ) 



k=l r 1=0 
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Rewriting the exponentials as cosines and taking limits gives the general identities. 



There is a slight difference in the form depending on the parity of p. For p even 



n 

k=i 1 



E 

1 < m < p — 1 
m odd 



2 cos ■ 



mx 



p n 



smx 



For p odd 



fc=i F 



i + 



E 



2 cos 



V 



1 < m < p — 1 
m even 



For p = 6 the identity takes the form 



nl / x ox ux 

« 2cos ^ + 2cos ^ + 2cos ^ 



fc=i 



3x 



5.x 



6 fc 



6 fc 



For p = 7 the identity takes the form 



n 

k=l 



2x Ax 6x 

1 + 2 cos + 2 cos + 2 cos 



smx 



smx 



smx 



(12) 



(13) 



(14) 



(15) 



For larger p fewer terms in the product are needed for the same degree of accuracy in the approxi- 
mation to sinx/x. In fact, by letting p go to infinity the first factor alone approaches sinx/x and 
provides a novel derivation of a well-known result. I leave it to the reader to work it out. 



4 Probabilistic Interpretation 



Mark Kac, in his delightful and now classic Carus monograph |2j, proves the first cosine identity 
Q in a way that is equivalent to the one we have outlined, although he does not explicitly use the 
Fourier transform, delta functions, and convolution. He then turns the identity into a question of 
probability, which for him was the leitmotif of his mathematical work. 

The original product identity Q arises from the following experiment. Flip a fair coin repeatedly. 
Beginning with 0, add 1/2 if the result is heads and subtract —1/2 if the result is tails. On the 
next toss add or subtract 1/4; on the next add or subtract 1/8, and so on. What is the distribution 
of the sums over the probability space whose elements are the countable sequences of coin tosses? 
Clearly, the sums are distributed uniformly between —1 and 1. 

Let s n denote the nth partial sum. It is a sum of independent random variables a\ + a<i + • • • + a n , 
where has the probability distribution (1/2) (S 1 / 2 k +<5_ 1 / 2 fc)- The probability distribution of 
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a sum of independent random variables is the convolution of the respective distributions of the 
random variables. Therefore, s n has the distribution 

n 1 1 

* n 2 ( 5i ' 2k + 5 - i ' 2h ) = 2» + " ■ + %0 • 

fc=l 



The inverse Fourier transform of a probability measure is called its characteristic function. Thus, 
the characteristic function for the distribution of s n is the product YYk=i cosx /^ k - ^ n the theory 
of probability and statistics, characteristic functions are a powerful tool. Typically computations 
are done with characteristic functions in order to draw conclusions about distributions of random 
variables as in the standard proof of the Central Limit Theorem. Here, however, we have inverted 
the relationship in order to compute with the probability measures and to get results about the 
characteristic functions. 



5 Related Products: Examples and Conjectures 



5.1 Coin tossing and Cantor sets 



The Cantor set K is the set of points between and 1 whose ternary expansion has no l's in it. So 
z is in K if z = Y^k=\ ^fc^ , ifc S {0,2}. Define K n to be the set of elements of K that have the 
form Ylk=i tk3~ k , arid define a probability measure supported on K n 



on ^ ] (16) 



z&K n 

K n has 2 n elements so jx n is equally distributed on K n . The sequence (/x n ) has a limit fj,, which 
can be described as assigning the following limit as the measure of a set E: 

KE) = hm t^JSlL. (17) 

n— >oo Z 

The measure fj, is also the Lebesgue-Stieltjes measure of the Cantor function. The Cantor function 
is continuous, non-decreasing, and has derivative zero on the complement of the Cantor set. Thus 
it defines a measure supported on the Cantor set, which is precisely the measure \i defined in (|17|). 

What is of interest in this note is that \i n is the finite convolution product 

f[^o + 5 2 /3*)- (18) 

k=l 

Consider the experiment of tossing a fair coin. On toss number k let 

_ heads 
Qk ~ i 2/3 fc tails 



k=l 
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Let s n = YHt=i a k- Then s n is equally distributed over K n . The characteristic function for the 
distribution of s n is IlLit 1 / 2 ) ( x + e 2 ^ 3 *). Define 

oo 

/w=n^( i + e2W3fc )- ^ 
fc=i 

(One checks easily that the product is convergent.) Then f = fx, the Cantor measure, but is it 
possible to characterize / in any other way? 

This leads us to look at the related infinite product Y\T=i cos Because 

•^(cos = i (<5 2/3 fc + <5_ 2/3fc ) 

the probabilistic interpretation is clear: add or subtract 2/3 fc on the kth toss with equal probability. 
Let s n be the sum of the first n values. What is the distribution of s n and what is the distribution 
of s = lim^oo s ra ? The exercise of expanding and plotting the values of S3 lead one to suspect that 
s is distributed "uniformly" over the Cantor set constructed from [—1,1] by successively removing 
middle thirds. That is easy to prove, as follows. 

Define the affine map of [0, 1] to [-1, 1] by z i-> 2(z- 1/2). Let z = E*fc3~ fc , t k £ {0,2}, be a point 
in the Cantor set. The ternary expansion of 1/2 is ^3 _fe , and so 2(z — 1/2) = ^ 2(tk — l)3 _fc . 
The coefficients 2(tk — 1) are either 2 or —2 with equal probability. 



This shows that the infinite product 11*^=1 cos §§ nas Fourier transform equal to the Cantor measure 
on the Cantor set constructed from [—1, 1] by removing middle thirds, but it does not give us a 
closed form like sinx/x. It would be most surprising if there were any simpler description of 
nfcLi cos §§- I n Figure 2 is a plot of the partial product with n = 8 and < x < 100. (The 
function is even.) Over this range the infinite product is indistinguishable from the eighth partial 
product. The self-similarity of the Cantor set at smaller and smaller scales appears to be reflected 
in the self-similarity of the graph at higher and higher frequencies. 




Figure 2. 
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5.2 Harmonic Series With Random Signs 



We have been looking at the sums of series of the form 

(20) 



k=l 

where t/~ is randomly chosen to be 1 or —1 with equal probability. Rademacher proved that if 
^2 Cfe 2 < co, then the sum converges with probability one on the probability space Q = {—1, 1} N . 
(CI can be identified with the unit interval and the probability measure with Lebesgue measure by 
using binary representations of numbers in the interval.) In j^j Kac gives the proof of this theorem 
due to Paley and Zygmund. It is also a theorem that the series diverges with probability one if 
^2 Cfc 2 = oo. Let us consider the random harmonic series 

Ep < 21 > 

k=l 

which converges almost surely by Rademacher's result, with the goal of understanding the distri- 
bution of the sums. This means we want to understand the distribution of the random variable 
s defined on $7. If we let s n be the partial sum, also a random variable, then the probability 
distribution of s n is the measure 

n 

»n = * II o (<Vfc + ( 22 ) 
k=l 



and its inverse transform is 



F- 1 (ri(x) = Hcos l . (23) 

k=l 

The product converges uniformly on compact sets as n — > oo, and so it is plausible that the sequence 
fi n converges to a probability measure [i that is the distribution of the random variable s. There 
is, however, a fair bit of analysis to make this rigorous. Assuming that the analysis can be made 
rigorous, then the plot of the Fourier transform of the infinite product IlfcLi cos § wm show how 
the sums are distributed. Let us call this function <P(uj). Then 

J -°° k=i 
i r°° °° x 

= — / (cos u)x + i sin ux) TT cos — dx (26) 
2tt l_ nn - LJ - k 

J 00 k=i 

1 f°° TT X 

= — coso;x cos — ax. (27) 

* Jo Li k 

There is not a closed form for 4>(u>) and so we resort to numerical integration. We truncated the 
infinite product at n = 1000 and integrated from to 15 using a straightforward Riemann sum 
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with dx = 0.02 and the midpoints of the subintervals for the points of evaluation. Values for to 
were from to 3.8 in multiples of 0.2. The integration was done with True BASIC on a portable 
Macintosh. See Figure 3. The distribution is very flat for — 1 < u < 1, much flatter than a normal 
distribution. A few of the computed values are given in this table 

] 





4>{uj) 


0.0 


.249995 


0.1 


.249991 


0.2 


.249972 


0.4 


.249809 


0.6 


.249092 


0.8 


.246819 


1.0 


.241289 


1.2 


.230494 


1.4 


.212941 


1.6 


.188425 


1.8 


.158271 


2.0 


.125000 


2.2 


.091729 


2.4 


.061576 


2.6 


.030596 


2.8 


.019506 


3.0 


.008711 


3.2 


.003181 


3.4 


.000908 


3.6 


.000192 


3.8 


.000028 



The value of ^>(0) is suspiciously close to 1/4, suggesting perhaps that 7r/4 is the value of the 
integral 



/ ]Tcos-cfe. (28) 
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One might also conjecture that Jq 00 cos 2x Y\T=i cos I ^ x = 7r /^- 

0.3 | , , , , , , 








-4-3-2-101234 



Figure 3. Graph of <f>{uS). 

For additional evidence we turned to simulations of the sums. Using MATLAB we ran 5000 sums 
of Y^k=i T with the values of tk picked randomly as ±1 with equal probability. Figure 4 shows a 
histogram of the sums. 

450 1 , , , , , , , 1 

400 - | 1 I — _ 

350 - — 1 

300 - 

250 - 

200 - 

150 - 

100 - 
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Figure 4. Histogram of 5000 random sums. 
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